Phases of Josephson Junction Ladders 



Colin Denniston 1 and Chao Tang 2 
1 Department of Physics, Princeton University, Princeton, New Jersey 08544 
2 NEC Research Institute, 4 Independence Way, Princeton, New Jersey 08540 

(February 1, 2008) 

Abstract 

We study a Josephson junction ladder in a magnetic field in the absence 

of charging effects via a transfer matrix formalism. The eigenvalues of the 

transfer matrix are found numerically, giving a determination of the different 

phases of the ladder. The spatial periodicity of the ground state exhibits 

a devil's staircase as a function of the magnetic flux filling factor /. If the 

transverse Josephson coupling is varied a continuous superconducting-normal 

transition in the transverse direction is observed, analogous to the breakdown 

of the KAM trajectories in dynamical systems. 
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Two-dimensional arrays of Josephson junctions have attracted much recent theoretical 
and experimental attention ]!]. Interesting physics arises as a result of competing vortex- 
vortex and vortex-lattice interactions. It is also considered to be a convenient experimental 
realization of the frustrated XY models. In this paper, we discuss the simplest such system, 
namely the Josephson junction ladder (JJL) [Q-^J shown in the inset of Fig. [TJ 

To construct the system, superconducting elements are placed at the ladder sites. Below 
the bulk superconducting-normal transition temperature, the state of each element is de- 
scribed by its charge and the phase of the superconducting wave function || . In this paper we 
neglect charging effects, which corresponds to the condition that 4e 2 / C <C J, with C being 
the capacitance of the element and J the Josephson coupling. Let 8j {6'-) denote the phase on 
the upper (lower) branch of the ladder at the j'th rung. The Hamiltonian for the array || can 
be written in terms the gauge invariant phase differences, jj = 9j — — (27r/0 o ) JJ_i A x dx, 
7 j = 0'. - Q'._ x - (2tt/0 o ) If -i A x dx, and a, = % - 6, - (2tt/0 o ) jf A y dx: 

H = - J2(J X cos jj + J x cos + J y cos atj), (1) 
j 

where A x and A y are the components of the magnetic vector potential along and transverse 
to the ladder, respectively, and 0o the flux quantum. The sum of the phase differences around 
a plaquette is constrained by jj — + ctj — ctj-i = 2n(f — rij), where rij = 0, ±1, ±2, ... 
is the vortex occupancy number and / = <f)/<f)o with <p being the magnetic flux through a 
plaquette. With this constraint, it is convenient to write Eq. ([[]) in the form 

Ti = — Jj^{2 cosr/j cos[(a,_i — aj)/2 + ir(f — rij)} 
j 

+ J t cos atj}, (2) 

where rjj = (7^ + 7j)/2, J = J x and J t = J y /J x - The Hamiltonian is symmetric under 
/—>•/ + 1 with rij — > rij + 1, and / — > — / with rij — > — rij, thus it is sufficient to 
study only the region < / < 0.5. Since in one dimension ordered phases occur only at zero 
temperature, the main interest is in the ground states of the ladder and the low temperature 
excitations. Note that in Eq. @ rjj decouples from <y,j and rij, so that all the ground states 
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have r/j = to minimize Ji. The ground states will be among the solutions to the current 
conservation equations dTC/daj = 0: 

Jt sinaj = sin[(otj_i — otj)/2 + ir(f — %)] 

- wa[(aj ~ «i+i)/2 + tt(/ - n j+1 )]. (3) 

For any given / there are a host of solutions to Eq. (H). The solution that minimizes the 
energy must be selected to obtain the ground state. 

If one expands the inter-plaquette coupling term in Eq. (§), cos[(<x,_i — aij)/2 + ir(f — rij)], 
about it's maximum, the discrete sine-Gordan model is obtained. A vortex (rij = 1) in the 
JJL corresponds to a kink in the sine-Gordan model. Kardar used this analogy to argue 
that this system should show similar behavior to the discrete sine-Gordan model which has 



been studied by several authors [§-10]. This analogy is only valid for J t very small so that 
the inter-plaquette term dominates the behavior of the system making the expansion about 
its maximum a reasonable assumption. However, much of the interesting behavior of the 
discrete sine-Gordan model occurs in regions of large Jt {Jt ~ 1)- Furthermore, much of 
the work by Aubry || on the sine-Gordan model relies on the convexity of the coupling 
potential which we do not have in the JJL. 

In this Letter we formulate the problem in terms of a transfer matrix obtained from 
the full partition function of the ladder. The eigenvalues and eigenfunctions of the transfer 
matrix are found numerically to determine the phases of the ladder as functions of / and J t . 
We find that the spatial periodicity of the ground states goes through a devil's staircase as a 
function of /. We then study the properties of various ground states and the low temperature 
excitations. As J t is varied, all incommensurate ground states undergo a superconducting- 
normal transition at certain J t which depends on /. One such transition will be analyzed. 
Finally we discuss the critical current. 

The partition function for the ladder, with periodic boundary conditions and K = 
J/ksT, is 
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/TV 
daidr)ie^p{K (2 cos cos[(aij_i - a i )/2 + n(f - n*)] + J t coso!j)} . (4) 

* {Ui} 

The ^ can be integrated out resulting in a simple transfer matrix formalism for the parti- 
tion function involving only the transverse phase differences: Z = Yu f-. K d a iP( a i-ii a i) = 
Tr P N . The transfer matrix elements P(a,a') are 

P(a,a') = 4irexp[KJ t (cosa + cosa')/2]I (2Kcos[(a-a')/2 + irf]), (5) 

where Jo is the zeroth order modified Bessel function. Note that the elements of P are real 
and positive, so that its largest eigenvalue Ao is real, positive and nondegenerate. However, 
since P is not symmetric (except for / = and f = 1/2) other eigenvalues can form complex 
conjugate pairs. As we will see from the correlation function, these complex eigenvalues 
determine the spatial periodicity of the ground states. 
The two point correlation function of ay's is 



( e <(ao-a,)) = lim 



(nf r^p(^-i,^))e^°-"') 



7V-*oo Z 
n V A 0/ 

where we have made use of the completeness of the left and right eigenfunctions. (Note 

that since P is not symmetric both right ip^ an d left ip^ eigenfunctions are need for the 

evaluation of correlation functions.) The A n in Eq. fl6|) are the eigenvalues (|A n | > |A n +i| 

and n — 0, 1,2, ...), and the constants c n = da'ip^ {a')e la ip^(a') dail)^{a)e~ ia ipQ{a). 

In the case where Ai is real and |Ai| > \X2\1 Eq. © simplifies for large I to 
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A<*o-*i) ) = Co + c (±L\ ; |A 1 |>|A 2 |. 



In the case where Ai = A 2 = |Ai|e* 27r= , Eq. (||) for large I is 
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(e l(ao " a,) ) = c + ( Cl e i27TEl + c 2 e~ i2nSl 



Ai 



, Ai — A2. 



I A c 

Note that while the correlation length is given by £ = [In | Ao/Ai|] _1 the quantity H = 
Arg(Xi)/27r determines the spatial periodicity of the state. Calculating A n numerically J?|, 



we found that for / smaller than a critical value f c \ which depends on J t , both Ai and 
A 2 are real. These two eigenvalues become degenerate at / c i, and then bifurcate into a 
complex conjugate pair. H as a function of / is shown in Fig. [I] for several different values 
of Jt- The shape of the curves in Fig. [I] is generally referred to as a devil's staircase. The 
steps of the staircase are at H = p/q, where p and q are integers. These are commensurate 
states with p vortices in each unit cell which consists of q plaquettes. For small Jt, the 
flat steps are connected by smooth differentiable curves; most states on the S — / curve 
are incommensurate states. As Jt increases, more and more steps appear and grow at the 
expense of the smooth regions. At J t = J f c ~ 0.7 the staircase becomes complete, i.e. there is 
a step for every rational H and the set of / which correspond to irrational H has zero measure. 
For J t > J^, the staircase becomes over-complete, i.e. steps of lower order rationals grow 
and steps of higher order rationals disappear ||11|| . Another important characterization of 
a state is the phase density p(a): p(a)da is the average fraction of all sites in the ladder 
with a < ai < a + da. If p(a) is a smooth function and p(a) > for a G (— 7T, ir] at 
T = 0, the ground state energy is invariant under an adiabatic change of a's. Consequently, 
there is no phase coherence between upper and lower branches of the ladder and hence no 
superconductivity in the transverse direction. In this case, we say that the a's are unpinned. 
If there exist finite intervals of a on which p(a) = 0, there will be phase coherence between 
the upper and lower branches and we say that the a's are pinned. In term of the transfer 



matrix, the phase density is the product of the left and right eigenfunctions of A |12 
p(a) = ip^(a)ip§(a). 

We first discuss the case where / < f c \. These are the "Meissner" states in the sense 
that there are no vortices (n, = 0) in the ladder. The ground state is simply ai = 0, 
7j = it f and 7^ = — 7vf, so that there is a global screening current ±J x simrf in the upper 
and lower branches of the ladder 0. The phase density p(a) = 5(a). Fig. |l](c) shows 
that at J t = 1, the Meissner state extends all the way from / = to / = f c \ ~ 0.28. 
The properties of the Meissner state can be studied by expanding Eq. (g) around = 0: 
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T~Cm = (^/4) J2j[ cos ( 7l f)( a j-i ~ a j) 2 + 2Jtaf]. The current conservation Eq. (|3|) becomes 



ttj+i = 2 (1 + J t / cos7r/) ctj — dj-\. 



(7) 



Besides the ground state ctj = 0, there are other two linearly independent solutions <x,- 
e ±i/£M Q f ^ which describe collective fluctuations about the ground state, where 



Cm 



1 + -^: + 



COS 7T j 



2J t 



f ~*~ \ COS7T/ ^COS7r/ 



Jt 



(8) 



£m is the low temperature correlation length for the Meissner state. (Note that Cm < 1 
for J t ~ 1 making a continuum approximation invalid.) As / increases, the Meissner state 
becomes unstable to the formation of vortices. A vortex is constructed by patching the two 
solutions of Eq. (|7|) together using a matching condition. The energy e v of a single vortex is 
found to be 



e v ss [2 + (tt 2 /8) tanh(l/2£ M )] costt/ 
- (tt + 1) sin 7t/ + 2J t , 



(9) 



for J t close to one. The zero of e v determines f c ± which is in good agreement with the 
numerical result from the transfer matrix. For f > f c i, e v is negative and vortices are spon- 
taneously created. When vortices are far apart their interaction is caused only by the expo- 
nentially small overlap. The corresponding repulsion energy is of the order Jexp(— 1/Cm)i 
where I is the distance between vortices. This leads to a free energy per plaquette of 
F = e v /l + Jexp(— /I 01- Minimizing this free energy as a function of / gives the 
vortex density for / > f c \\ (rij) = l" 1 = [£jtfhi |/ci — where a linear approximation is 
used for / close to f c \. 

We now discuss the commensurate vortex states, taking the one with H = 1/2 as an 
example. This state has many similarities to the Meissner state but some important differ- 
ences. The ground state is 
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a = arctan 



Jt 



sin(vr/) 



Oil = -"0, Oi i±2 = Oif, 



n Q = 0, m = 1, n i±2 = Hi, 



(10) 



so that there is a global screening current in the upper and lower branches of the ladder 
of ±2nJ(f — l/2)/y / 4 + J?. The existance of the global screening, which is absent in an 
infinite 2D array, is the key reason for the existance of the steps at S = p/q. It is easy to see 
that the symmetry of this vortex state is that of the (antiferromagnetic) Ising model. The 
ground state is two-fold degenerate. The low temperature excitations are domain boundaries 
between the two degenerate ground states. The energy of the domain boundary Jet, can be 
estimated using similar methods to those used to derive Eq. @ for the Meissner state. 
We found that e& = e° — (tc 2 /^A + Ji)\f — 1/2 1 , where e° depends only on J t . Thus the 
correlation length diverges with temperature as £ ~ exp(2 Je^/Zc^T). The transition from 
the S = 1/2 state to nearby vortex states happens when / is such that e& = 0; it is similar to 
the transition from the Meissner state to its nearby vortex states. All other steps H = p/q 
can be analyzed similarly. For comparison, we have evaluated £ for various values of / and 
T from the transfer matrix and found that £ fits £ ~ exp(2Jeb//csT) (typically over several 
decades) at low temperature. The value of e b as a function of / is shown in Fig. ^|for J t = 1. 
The agreement with the above estimate for the H = 1/2 step is excellent. The tips of the 
peaks in Fig. |2] for states with H = \/q fit the relationship r ~ exp(— q/lo) with l ps 0.77, 
which is in good agreement with £m(/ = 0) ~ 0.76 of Eq. (^j). 

We now discuss the superconducting-normal transition in the transverse direction. For 
J t = 0, the ground state has 7$ = 7^ = and 

«i = 27r/i + «o -27T Yl'^o n i- ( U ) 

The average vortex density (rij) is /; there is no screening of the magnetic field. «o in 
Eq. (|Tl| ) is arbitrary; the ct's are unpinned for all /. The system is simply two uncoupled 
ID XY chains, so that the correlation length £ = 1/ksT. The system is superconducting 
at zero temperature along the ladder, but not in the transverse direction. As J t rises above 
zero we observe a distinct difference between the system at rational and irrational values of 
/. For / rational, the a's become pinned for J t > (p(a) is a finite sum of delta functions) 
and the ladder is superconducting in both the longitudinal and transverse directions at zero 
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temperature. The behavior for irrational / is illustrated in the following for the state with 
H = a g , where a g ~ 0.381966 ■ ■ • is one minus the inverse of the golden mean. Fig. ^ displays 
p(a) for several different Jt at H = a g . We see that the zero- frequency phonon mode (the 
smoothness of p(a)) persists for small J t > until a critical value Jt(f) ~ 0.7 where the a's 
become pinned and the ladder becomes superconducting in the transverse direction. In the 
sine-Gordan model, the pinning transition of this state coincides with the devil's staircase of 
Fig. |l| becoming complete |§J§ (If the Oj's are pinned in this state, then all incommensurate 
states should be pinned). The pinning transition of the incommensurate states can be also 
studied using Eq. (|3|) which can be viewed as a two-dimensional map. The disappearance of 
the zero-frequency phonon mode for irrational S's at finite small Jf{f) is equivalent to the 
breakdown of the Kolmogorov- Arnold- Moser (KAM) trajectories of the map ||13|| . 

We now turn to the subject of critical currents along the ladder. One can obtain an 
estimate for the critical current by performing a perturbation expansion around the ground 
state (i.e. {rij} remain fixed) and imposing the current constraint of sin 7.,- +sin7^ = /. Let 
8jj, Sjj and 5a j be the change of jj, 7'. and ctj in the current carrying state. One finds that 
stability of the ground state requires that Saj = 0, and consequently 8jj = Sjj = 1/2 cos jj. 
The critical current can be estimated by the requirement that the jj do not pass through 
7r/2, which gives J c = 2(7r/2 — 7 max ) cos7 max , where 7 max = max, (77). In all ground states 
we examined, commensurate and incommensurate, we found that 7 max < 7r /2, implying a 
finite critical current for all /. 

In conclusion, we have studied the equilibrium behavior of a Josephson junction ladder in 
a magnetic field in the absence of charging effects. The screening current plays an important 
role in this system. For f < f c i, there is a "Meissner state" with no vortices. For f > f&, 
the spatial periodicity of the ground state climbs a devil's staircase as a function of /. 
All incommesurate states undergo a superconducting-normal transition in the transverse 
direction as J t is increased, so that for J t > J t c ~ 0.7 the ladder is superconducting in 
both the longitudinal and transverse directions for all /. The critical current along the 
ladder is found to be finite for all /. Finally, although in one dimension there is no finite 

8 



temperature phase transition and no true long range order, our study showed that the 
correlation lengths in vortex states are extremely long for reasonably low temperatures. 
Thus it would be interesting to test these results experimentally. More extensive details 
of this system, including current carrying states will be presented elsewhere along with a 
similar analysis of the sine-Gordan model J7|. 

We thank Sue Coppersmith, Xinsheng Ling, and Qian Niu for many valuable discussions. 
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FIGURES 

FIG. 1. Periodicity, H = Ar 5 (Ai)/27r versus / for k B T/J = 0.005 and, (a) J t = 0.3, (b) J t 
= 0.7, (c) Jt = 1.0. Inset: The Josephson junction ladder is formed by the arrangement of the 
superconducting islands. The field H is out of the page and the arrows indicate the direction of 
the gauge invariant phase differences. 

FIG. 2. Effective Ising coupling as a function of / for J t = 1. The inset shows the statistical 
error for in the fitting. 

FIG. 3. p(a) = i/j^(a)tp^(a) versus a at k B T/J = 0.004 and E = 0.381966 011 • • •, and for (a) 
Jt = 0.4; (b) Jt = 0.65; (c) J t = 0.7; and (d) J t = 0.9. Note the smaller scale for the upper plots. 
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